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This paper addressses the kinetics and dynamics of a fam- 
ily of domain wall solutions along classical antiferromagnetic 
Heisenberg spin chains at low energies. The equation of mo- 
tion is derived and found to have long range position- and 
velocity-dependent two-body forces. A 'quiescent' regime is 
identified where the forces between walls are all repelling. 
Outside this regime some of the interactions are attractive, 
giving rise to wall collisions whereupon the colliding walls 
annihilate. The momentum of the system is found to be con- 
served in the quiescent regime and to suffer discontinuous 
jumps upon annihilation. The dynamics are illustrated by 
an exact solution for a double wall system and a numerical 
solution for a many-wall system. 

On circular chains the equations support stable periodic 
stripes that can rotate as a rigid body. It is found that the 
stripes are more stable the faster they rotate. The periodic 
structure can be destabilised by perturbing the walls' angular- 
velocity in which case there is a transition to another periodic 
structure, possibly via a cascade of annihilation events. 

75.60.Ch 



Spin chains with antiferromagnetic Heisenberg ex- 
change coupling have attracted much attention in recent 
years due to the remarkable progress in synthesis meth- 
ods, new experimental techniques, and a variety of pos- 
sible applications. These chains exhibit a rich behaviour, 
much of which is not fully understood, either in the quan- 
tum or the classical regimes Jj] . Of particular interest is 
the low-energy regime which is amenable to some analyt- 
ical treatment. It has been known for some time that this 
regime accommodates spin waves and instanton solutions 
H H . Recently, however, it has been shown to also sup- 
port a family of multitwist solutions Q] || that extend a 
single twist found earlier || . These solutions represent an 
arbitrarily number of domain walls that can travel along 
the chain. The multitwists are found by starting from the 
discrete Hamiltonian for the exchange coupling between 
nearest neighbour spins on the antiferromagnetic chain 
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and distinguishing between the even, S e , and odd, S , 
sublattices. Defining the total, ( oc S e + S , and stag- 
gered, ff oc S e — S magnetisations (properly normalised) 
and taking the continuum limit, leads in the low energy 
regime to 
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where t is renormalised time and x is the arclength pa- 
rameter along the chain. This equation translates into 
two coupled equations for the polar and azimuthal an- 
gles, 9 and <j>, of the unit vector ff which, in turn, are 
identical to the so called Belavin-Polyakov equations || 
in the context of the nonlinear sigma model. The resem- 
blance is not complete because in this context the equa- 
tions are derived via the Bogomol'nyi derivation and 
consequently come complete with uniform boundary con- 
ditions. Not do in the context of the antiferromagnetic 
chain where the boundary conditions can be arbitrary. 

Introducing cos 9 = tanh0, whereupon ff = 
(cos<^>, sin^, sinh^))/cosh^> Q, the coupled equations are 
transformed into 

dtp def) dip def) 
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In the absence of sources, which give rise to logarithmic 
terms (corresponding to n-instanton solutions || ||), the 
solutions are all the analytic functions in the x — t plane. 
The smooth behaviour of these functions, once operated 
on with the hyperbolic function, yields multitwists in the 
angle 9 along the chain where the nodes of ip are. These 
multitwists correspond to domain walls in the staggered 
magnetisation ff. 

The aim of this paper is to analyse the dynamics of 
the domain walls corresponding to these twists. Specif- 
ically, I will derive the equation of motion of the walls, 
showing that they interact via long range terms that are 
position and velocity dependent. It will be argued that 
there is a 'quiescent' regime where the system of walls 
conserves its global momentum. When wall velocities 
become sufficiently high the system exits this regime and 
domain walls develop mutual attarction, leading to po- 
tential annihilation events and consequently to rich dy- 
namics. Upon annihilation, momentum conservation is 
broken for the entire system but remains intact for the 
remaining subsystem. These results will be demonstrated 
on a double- and a many- wall systems. Next, I will show 
that closing the chain to form a ring gives rise to a fam- 
ily of periodic stripe solutions in the quiescent regime. 
These lattice structures are stable to small perturbations 
but large perturbation effect annihilations, which lead to 
transitions to other periodic solutions. It will be also 
argued that a cascade of such transitions may occur re- 
sulting in intriguing dynamics. 

Consider an open-ended infinite chain prepared with 
an initial set of domain walls that occupy a finite region. 
The general analytic solution for ip can be represented in 
the form 
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Where N is the total number of domain walls determined 
by the initial data. The coefficients P m are polynomials 
of order m in t whose structure has been discussed in ref- 
erence 0]. Real roots, x n correspond to physical domain 
walls in the field if where the angle 9 goes through an 
abrupt change of tt. The aim in the following is to find 
and analyse the dynamics of these walls. 

First note that the field -0 must satisfy Laplace's equa- 
tion and therefore, carrying out the explicit derivatives, 
we get 
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where x n is the location of the nth wall. The solution for 
if) is valid for all x and in particular when x approaches 
the location of the nth wall, x n . It follows that the 
residues of the poles at x n must all vanish. A straight- 
forward calculation yields the equation of motion for the 
nth wall 
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This innocent-looking set of equations turns out to de- 
scribe very rich dynamics and the rest of the paper is 
devoted to an analysis of its physical consequences. 

There are several features to be noted about the set 
of equations @: (1) The force on a wall is the sum of 
two-body forces between it and all the other walls along 
the chain. These forces are long ranged and depend on 
both position and velocities. (2) Denoting the force on 



wall n due to wall m by f nm then /„ 



-fmn, which is 



reminiscent of Newton's third law of action and reaction 
in translational invariant systems that this one is not. (3) 
As a consequence of the above property, the system of 
walls conserves momentum. To see this define the centre 
of mass, 7 = 2fcLi X k/N and summ over n. The terms 
in the sum over m cancel in pairs leading to 7 = and 
hence to 7 = po — constant. (4) Eqs (H) are symmetric 
under time reversal tt — > —t. Therefore, for each forward 
solution ij}f(t) there exists a bacward one ipb(t) — i/)f(—t). 

To illustrate the kinetics consider a double wall system. 
The equations of motion read 
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Defining the centre of mass velocity 7 = (xi + x 2 ) /2 = pa 
and their separation S — (x 2 — xi)/2, one arrives at an 
equation for 8 
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where c = \J(1 + pf, ) ■ This equation is symmetric not 
only under time reversal, but also under i5 — > —8. A 
particular solution of this equation is 
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where the positive (negative) sign corresponds to a so- 
lution that propagates forward (backward) in time. To 
find the general solution define 8 = c rj 2 / 3, and substitute 
into eq @, which gives 
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Integrating this equation we have 
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with E the integration constant. This equation can now 
be solved by quadratures 
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where \ — AcE/9. Again the positive (negative) sign 
corresponds to forward (backward) propagation, which 
increases (decreases) monotonically with time. The case 
E = yields the linear solution (|l0|). There can be three 
scenarios for the motion of the two walls: (a) moving 
apart, (b) approaching each other and eventually col- 
liding, or (c) first approaching and then moving away. 
Which scenario takes place depends sensitively on initial 
conditions. Solutions (|l^) and ( |l3| ) are real for all E > 
when 8 > and represent either the two walls moving 
apart forever (forward propagation) or moving into an 
eventual collision (backward propagation). In the latter 
case the time to collision is t co u — 8q/c for the linear 
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eral case. A careful analysis of E < gives that in this 
case the walls first approach each other until delta = \ X \ 
whereupon the solution switches to another branch and 
they move back out. This happens when the approach 
is at sufficiently low velocities such that their interaction 
is always repelling. Once the walls have collided in the 
second scenario they annihilate. Whether they will reap- 
pear at a later time can be figured out from eq (0) for 
ip. It is worth noticing that on forward propagation the 
separation between the walls converges asymptotically to 
a steady increase at a rate of 8 = c that is independent of 
all initial conditions except the total momentum of the 
system pq . 

Having understood the double wall case let us turn to 
many-wall systems. Consider two walls in such a sys- 
tem, say n and to, whose velocities are, respectively, x n 
and x m - The interaction between these walls is repul- 
sive when Xnim > — 1 and is attractive otherwise. The 
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FIG. 1. The dynamics of ten walls along the antiferro- 
magnetic spin chain in the quiescent regime. The walls initial 
positions are normally distributed around an arbitrary point 
along of the chain and they are assigned low velocities chosen 
from a Gaussian distribution of zero mean. Note how the ini- 
tial closure of the gap between some walls is reversed due to 
walls repulsion. 

attraction comes about when the two walls either move 
away from, or towards, each other at sufficiently high 
speeds. In the former case the attraction slows down too 
fast a separation, as is demonstrated in the double wall 
system where 5 — > c. In the latter case, attraction be- 
tween walls that approach each other too fast is bound 
to lead to a collision unless deflected by other walls in 
the system. On slow approaches walls do not collide but 
rather achieve a minimal separation before departing in 
opposite directions. The regime where all pairs satisfy 
i n x m > —1 is 'quiescent' in the sense that no walls can 
collide and annihilate. The number of walls in this regime 
remains constant and therefore so does the total momen- 
tum. Once the system exits the quiescent regime the 
conservation laws do not hold. Upon each annihilation 
event the number of walls changes and the system's mo- 
mentum undergoes a discontinuous jump. The momen- 
tum lost in such a jump is exactly that of the two anni- 
hilating walls and it is the momentum of the remaining 
subsystem which remains conserved from the moment of 
annihilation onward. 

This analysis is supported by numerical studies where 
the above kinematic scenarios can be clearly observed in 
the dynamics of the walls. Figure 1 shows a system of ten 
domain walls in the quiescent regime. The walls initial 
positions are randomly chosen from a Gaussian distribu- 
tion around an arbitrary location along the chain. Their 
initial velocities are chosen from a Gaussian distribution 
of zero mean and a standard deviation of 0.5 (arbitrary 
units). As expected, the walls fan out with time, demon- 
strating the repulsive interactions. It can be observed 
that some neighbouring walls initially start approaching 
each other but their momentum is too low to effect colli- 
sion and eventually they pull apart, as discussed above. 
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FIG. 2. The dynamics of ten walls along the antiferro- 
magnetic spin chain outside the quiescent regime. The walls 
are initially positioned as in figure 1, but their random ini- 
tial velocities are chosen from a Gaussian distribution that is 
three times as wide. Note the collision and annihilation of 
two neighbouring walls. Before (after) annihilation the total 
momentum of the system is constant at -0.39 (-1.72) with the 
change occuring abruptly at the moment of annihilation. 

Figure 2 shows those walls starting from the same initial 
positions but at initial velocities that are chosen from a 
Gaussian distribution whose standard deviation is three 
times as large. In this case some walls (eg, the third 
and fourth from the right) start off approaching at ve- 
locities that are sufficiently high to give rise to collision, 
upon which they collide and annihilate. Computing the 
momentum of the system indeed gives that these are con- 
stant before and after annihilation at respective values of 
-0.39 and -1.72 (arbitrary units) with the change occuring 
abruptly at the moment of annihilation. 

Let us turn now to finite chains that are closed to form 
rings. Examination of the equations of motion shows that 
they support periodic solutions of N walls around the 
ring that are located at intervals of a = L/N, where L is 
the chain length. These are stripe solutions that are sta- 
ble in the quiescent regime. This can be seen by applying 
a perturbation to one of the walls in this state, moving 
it from its lattice position by a distance smaller than a 
while keeping its angular velocity low enough so as not 
to exit the quiescent regime. The repulsive interaction 
with the closer neighbour dominates and pushes the wall 
back to its lattice position, as can be readily verified by 
linear perturbation analysis. This means that in the qui- 
escent regime the stripe lattice dynamics consists only of 
oscillations around the periodic structure. The spectrum 
of these oscillations would be interesting to analyse but 
this is outside the thrust of this paper. 

A most intriguing feature of these solutions is that ro- 
tating the entire stripe lattice rigidly at an arbitrary an- 
gular velocity u> is also a solution. Moreover, the ampli- 
tude of the repulsive forces between the walls increases 
with the angular velocity as 1 + (ojL/2tt) 2 . The rotat- 
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ing solution is again stable under a perturbation of the 
walls from their lattice positions. However, the stability 
of the periodic structure depends more on perturbations 
to their angular velocities then positions. Suppose we 
give wall n an excessive angular velocity Suj. If Suj is 
in the same direction as uj this wall experiences a net 
restoring force from its two neighbours. The same ap- 
plies if Suj is in the opposite direction to u> but its mag- 



nitude is smaller than 



Being stable, 
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the dynamics of the lattice in the rotating frame due to 
low kinetic energy excitations will again consist of oscil- 
lations of the walls around the lattice positions. Once 5u> 
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in the 



exceeds the critical value of uj 

direction, the net force on the nth wall becomes attrac- 
tive and pulls it further in that direction. The nth wall 
then runs into the (n — l)th wall and they annihilate. 
Thus, too high excessive angular velocities to individual 
walls destabilises the lattice. 

The dynamics ensuing in the wake of such an insta- 
bility are fascinating: Upon annihilation a gap appears 
in the periodic structure and the remaining N — 2 walls 
must shift to accomodate the new conditions. If this 
readjustment gives rise only to low excessive angular ve- 
locities then the system settles, without further ado, into 
a new structure of period a = L/(N — 2). If, however, 
the readjustment gives rise to high excessive velocities in 
the — uj direction then another annihilation may occur, 
leading to further damage of the lattice. Thus, a cascade 
of annihilation events may follow the initial instability 
leading to strange and rich dynamics. This author be- 
lieves that such behaviour can be observed by monitor- 
ing directly the spectrum of the stripe lattice to vibra- 
tional excitations. At low kinetic energies || one expects 
plasmon-like excitations due to the long range Coulomb 
interaction. But, at higher kinetic energies the spectrum 
should become muddled by the above dynamics. 

To conclude, the dynamics of domain walls on classi- 
cal antiferromagnetic Heisenberg chain have been found 
to be surprisingly rich. The equation of motion of the 
walls along the chain has been derived from first prin- 
ciples, showing that the walls experience only two-body 
forces. These forces consist of both a Coulomb-like term 
and a long range velocity dependent term. The interac- 
tions can be either repulsive, in a quiescent regime, or 
attractive at high velocities. The system of walls has 
been found to conserve momentum as long as the num- 
ber of walls is conserved, which is always the case in the 
quiescent regime. Outside this regime walls may collide 
and annihilate, upon which the total momentum of the 
system changes discontinuously to a new level that re- 
main constant until another annihilation event occurs. 
An explicit solution of a double wall system has been 
presented, giving insight into the dynamics of many-wall 
systems. A numerical solution of the latter has been 
discussed both inside and outside the quiescent regime 
and annihilation events have been demonstrated with in- 



creasing wall velocities. This solution also demonstrated 
the above principle of momentum conservation between 
annihilations. 

On circular chains the equations of motion give rise 
to periodic stripe solutions that are stable for small per- 
turbations. These periodic structures can roitate as a 
rigid body with the structures becoming more stable with 
increasing angular velocity. This observation may have 
direct application in rotating magnetic memory devices. 
The dynamics of the stripe lattice has been found to be 
very rich. In the quiescent regime it oscillates and at low 
kinetic energies the system should resemble a confined 
electron gas due to the long range Coulomb-like repul- 
sion. Increasing the kinetc energy brings into play the 
nonlinearity of the equation of motion distorting the os- 
cillation spectrum into potentially new and interesting 
configurations. Nevertheless, the dynamics should still 
be confined to oscilatory vibrational excitations. As the 
kinetc energies increase further the system exits the qui- 
escent regime, interactions turn attractive and walls can 
collide and annihilate. Annihilations on a ring lead to re- 
arrangements of the walls, which in turn, can either settle 
into a new periodic structure or result in a cascade of fur- 
ther annihilation events, the dynamics of which should be 
very interesting to study. 
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